Evidence for nonhadronic degrees of freedom in the transverse mass spectra of kaons 

from relativistic nucleus-nucleus collisions? 
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We investigate transverse hadron spectra from relativistic nucleus-nucleus collisions which reflect 
important aspects of the dynamics - such as the generation of pressure - in the hot and dense 
zone formed in the early phase of the reaction. Our analysis is performed within two independent 
transport approaches (HSD and UrQMD) that are based on quark, diquark, string and hadronic 
degrees of freedom. Both transport models show their reliability for elementary pp as well as light-ion 
(C+C, Si+Si) reactions. However, for central Au+Au (Pb+Pb) collisions at bombarding energies 
above ~ 5 A- GeV the measured K transverse mass spectra have a larger inverse slope parameter 
than expected from the calculation. Thus the pressure generated by hadronic interactions in the 
transport models above ~ 5 A- GeV is lower than observed in the experimental data. This finding 
shows that the additional pressure - as expected from lattice QCD calculations at finite quark 
chemical potential and temperature - is generated by strong partonic interactions in the early phase 
of central Au+Au (Pb+Pb) collisions. 

PACS numbers: 25.75.-q, 25.75.Dw, 25.75.Ld, 13.60.Le 



The phase transition from partonic degrees of freedom 
(quarks and gluons) to interacting hadrons is a central 
topic of modern high-energy physics. In order to un- 
derstand the dynamics and relevant scales of this transi- 
tion laboratory experiments under controlled conditions 
are presently performed with ultra-relativistic nucleus- 
nucleus collisions. Hadronic spectra and relative hadron 
abundancies from these experiments reflect important as- 
pects of the dynamics in the hot and dense zone formed in 
the early phase of the reaction. Furthermore, as has been 
proposed early by Rafelski and Miiller the strangeness 
degree of freedom might play an important role in dis- 
tinguishing hadronic and partonic dynamics. 

Estimates based on the Bjorken formula |2| for the en- 
ergy density achieved in central Au+Au collisions sug- 
gest that the critical energy density for the formation of 
a quark-gluon plasma (QGP) is by far exceeded during a 
few fm/c in the initial phase of the collision at Relativis- 
tic Heavy Ion Collider (RHIC) energies @, but sufficient 
energy densities (~ 0.7-1 GeV/fm 3 jj]) might already be 
achieved at AlternatmgGradient Synchrotron (AGS) en- 
ergies of ~ 10 A-GeV yjla]- More recently, lattice QCD 
calculations at finite temperature and quark chemical po- 
tential n q show a rapid increase of the thermodynamic 
pressure P with temperature above the critical temper- 
ature T c for a phase transition to the QGP. The crucial 
question is, however, at what bombarding energies the 
conditions for the phase transition arc fulfilled. 

Presently, transverse mass (or momentum) spectra of 
hadrons are in the center of interest: the measured trans- 
verse mass hit spectra of kaons 
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at AGS and SPS (from 30 to 160 A-GeV at the Super Pro- 
ton Synchrotron) energies show a substantial flattening 
or hardening of the spectra in central Au+Au collisions 
relative to pp interactions (cf. d, @)- This hardening 
of the spectra (or increase of the inverse slope parame- 
ter T in ||TJ) is commonly attributed to strong collective 
flow which is absent in the respective pp data. We will 
demonstrate in this Letter that the pressure needed to 
generate such a large collective flow cannot be produced 
by the interactions of hadrons in the expansion phase of 
the hadronic fireball. 

Two independent transport models that employ 
hadronic and string degrees of freedom, i.e. UrQMD 
El 113 and HSD [H El, are used for our study. They 
take into account the formation and multiple rescatter- 
ing of hadrons and thus dynamically describe the gener- 
ation of pressure in the hadronic expansion phase. The 
UrQMD transport approach El EI includes all baryonic 
resonances up to masses of 2 GeV as well as mesonic 
resonances up to 1.9 GeV as tabulated by the Particle 
Data Group Q. For hadronic continuum excitations 
a string model is used with hadron formation times in 
the order of 1-2 fm/c depending on the momentum and 
energy of the created hadron. The transport approach 
is matched to reproduce the nucleon-nucleon, meson- 
nucleon and meson-meson cross section data in a wide 
kinematic range |lCt 1 1 11] - In the HSD approach nuclc- 
ons, A's, N*(1440), N*(1535), A, £ and £* hyperons, 
S's, S*'s and f2's as well as their antiparticles are in- 
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FIG. 1: Comparison of transverse mass spectra from HSD 
(solid lines) and UrQMD (dashed lines) for tt~,K + and 
TW(xO.l) from 17.5% central C+C (upper part) and 12.5% 
central Si+Si (lower part) reactions at 160 A-GeV close to 
midrapidity with data from Ref. [Tfll ] (full symbols). 



eluded on the baryonic side whereas the CP and 1~ octet 
states are included in the mesonic sector. High energy 
inelastic hadron-hadron collisions in HSD are described 
by the FRITIOF string model 0] whereas low energy 
hadron-hadron collisions are modeled based on experi- 
mental cross sections. We point out importantly, that 
no parton-parton scattering processes are included in 
the studies below contrary to the multi-phase transport 
model (AMPT) |l6j |. which is currently employed from 
upper SPS to RHIC energies. 

Whereas the underlying concepts of the two transport 
theoretical models are very similar, the actual realiza- 
tions differ considerably Indeed, a systematic analysis 
of results from both models and experimental data for 
central nucleus- nucleus collisions from 2 to 160 AGeV 
has shown that the 'longitudinal' rapidity distributions of 
protons, pions, kaons, antikaons and hyperons are quite 
similar in both models and in reasonable agreement with 
available data. The exception are the pion rapidity spec- 
tra at the highest AGS energy and lower SPS energies, 
which are overestimated by both models ^3 • For a com- 
parison of HSD calculations with experimental data at 
RHIC energies we refer the reader to Ref. |Tgj . 

We focus on transverse mass spectra of pions and 
kaons/antikaons from central Au+Au (Pb+Pb) collisions 
from 2 A-GeV to 21.3 ATeV and compare to recent data. 
We start with results on nix spectra for central C+C and 
Si+Si collisions at 160 A-GeV. In Fig. 2] the experimen- 
tal ttit spectra for 7r~, K + and K~ for C+C reactions 
(17.5% centrality; upper part) and Si+Si (12.5% central- 
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FIG. 2: Comparison of transverse mass spectra from HSD 
(solid lines) and UrQMD (dashed lines) for -k^ , K + and 
K~(x0.1) from central Au+Au (Pb+Pb) reactions at 4, 30, 
and 160 A-GeV at midrapidity with the data from Refs. H,|2(| 
(full symbols). 

ity; lower part) 0] are compared to HSD (solid lines) 
and UrQMD calculations (dashed lines). The mr spec- 
tra are calculated in a rapidity interval (3.0 < yi a b < 
3.23) close to midrapidity. The agreement between the 
transport calculations and the data is quite satisfactory; 
no obvious traces of 'new' physics are visible. The situ- 
ation, however, changes for central Au+Au (or Pb+Pb) 
collisions as demonstrated in Fig. [3 where we display 
the calculated spectra from UrQMD (dashed lines) and 
HSD (solid lines) for 5% (or 7 %) central Au+Au (or 
Pb+Pb) collisions at 4, 30 and 160 A-GeV in compari- 
son to midrapidity data from Refs. |8|,|2Cj 32]. Whereas 
at the lowest energy of 4 A-GeV the agreement between 
the transport approaches and the data is still acceptable, 
severe deviations are visible in the if ± spectra at SPS 
energies of 30 and 160 A-GeV. Note that the ^ spectra 
are reasonably described at all energies while the inverse 
slope T of the K ± transverse mass spectra is underesti- 
mated severely by about the same amount in both trans- 
port approaches (within statistics). The increase of the 
inverse slopes in heavy-ion collisions with respect 
to pp collisions, which is generated by rescatterings of 
produced hadrons in the transport models, is only small 
because the elastic meson-baryon scattering is strongly 
forward peaked and therefore gives little additional trans- 
verse momentum at midrapidity. 

The question remains whether the discrepancies shown 
in Fig. |21 might be due to conventional hadronic medium 
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effects. In fact, the tot slopes of kaons and antikaons 
at SIS energies (1.5 to 2 A- GeV) were found to differ 
significantly j2l|. As argued in the different slopes 
could be traced back to repulsive potentials of kaons with 
nucleons, which lead to a hardening of the K + spectra, 
and attractive antikaon-nucleon potentials, which lead to 
a softening of the K~ spectra. However, the effect of such 
potentials was calculated within HSD and was found to 
be of minor importance at AGS and SPS energies 
since the meson densities are comparable to or even larger 
than the baryon densities at AGS energies and above. 

Additional self energy contributions stem from in- 
teractions with mesons; however, s-wave kaon-pion inter- 
actions are weak due to chiral symmetry arguments and 
p-wave interactions such as n + K «-> K* transitions are 
suppressed substantially by the approximately 'thermal' 
pion spectrum. A recent study on the kaon potentials in 
hot pion matter gives kaon mass shifts of about —52 MeV 
and vector potentials of ~ +49 MeV .22] for a pion gas 
at temperature T= 170 MeV. We have employed even 
slightly larger K ± potentials in dynamical HSD calcula- 
tions and achieved a hardening of the spectra of less 
than 10%. 

Furthermore, we have pursued the idea of Refs. [23ll24| 
that the spectra could be hardened by string-string 
interactions, which increase the effective string tension a 
and thus the probability to produce mesons at high m<r 
[Icil^ . In order to estimate the largest possible effect of 
string-string interactions we have assumed that for two 
overlapping strings the string tension a is increased by 
a factor of two, for three overlapping strings by a factor 
of three etc. Here the overlap of strings is defined geo- 
metrically assuming a transverse string radius R s , which 
according to the studies in Ref. [23 should be R s < 0.25 
fm. Based on these assumptions (and i? s =0.25 fm), we 
find only a small increase of the inverse slope parameters 
at AGS energies, where the string densities are low. At 
160 A- GeV, the model gives a significant hardening of 
the spectra by about 15%, which, however, is still much 
less than in the data in Fig. [21 

Our findings are summarized in Fig. [3| where the de- 
pendence of the inverse slope parameter T (see Eq. 
on \f s is shown and compared to experimental data 
[il Ha] for central Au+Au (Pb+Pb) collisions (l.h.s.) and 
pp reactions (r.h.s.). The upper and lower solid lines 
(with open circles) on the l.h.s. in Fig. [3 correspond to 
results from HSD calculations, where the upper and lower 
limits are due to fitting the slope T itself, an uncertainty 
in the repulsive i^-pion potential or the possible effect 
of string overlaps. The open triangles represent slope pa- 
rameters from default UrQMD calculations (version 1.3) 
which are well within the limits obtained from the dif- 
ferent HSD calculations. It's worth to point out that in 
the pure cascade mode the kaon slope T is slightly larger 
in the UrQMD calculations than from HSD at AGS and 
SPS energies, which is related to the isotropic decay of 
higher mass (< 2 GeV) resonances in UrQMD as well 
as a different string fragmentation function. The slope 
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FIG. 3: Comparison of the inverse slope parameters for K + 
and K~ mesons from central Au+Au (Pb+Pb) collisions 
(l.h.s.) and pp reactions (r.h.s.) as a function of the invari- 
ant energy y^s from HSD (upper and lower solid lines) and 
UrQMD (open triangles) with data from Refs. |a. |19|. The 
upper and lower solid lines result from different limits of the 
HSD calculations as discussed in the text. 



parameters from pp collisions (r.h.s. in Fig. [3J are seen 
to increase smoothly with energy both in the experiment 
(full squares) and in the HSD calculations (full lines with 
open circles). The error bars related to our calculations 
are due to an uncertainty in extracting the slope param- 
eter when considering different m-r intervals. 

We mention that the RQMD model |2^| gives higher 
inverse slope parameters for kaons at AGS and SPS en- 
ergies than HSD and UrQMD, which essentially can be 
traced back to the implementation of resonances with 
masses above 2 GeV as well as 'color ropes' that decay 
isotropically in their rest frame 26]. However, we con- 
sider the implementation of such 'heavy resonances' with 
unknown properties as a phenomenological fit to experi- 
mental data. 

This still leaves us with the question of the origin of 
the rapid increase of the slopes with invariant energy 
for central Au+Au collisions, which is missed in both 
transport approaches. We recall that higher transverse 
particle momenta either arise from repulsive self ener- 
gies in mean-field dynamics - or from collisions, which 
reduce longitudinal momenta in favor of transverse mo- 
menta. As shown in Fig. 3 conventional hadron self- 
energy effects and hadronic binary collisions are insuffi- 
cient to describe the dramatic increase of the slopes 
as a function of yfs. We conclude that this points to- 
wards a different mechanism for the generation of the 
pressure observed experimentally and propose, that par- 
tonic degrees of freedom should be responsible for this 
effect already at ~ 5 A-GeV. 

Our arguments are based on a comparison of the trans- 
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verse pressure from the transport models with recent 
lattice QCD calculations 0: The transverse pressure 
(P xx w Pyy) - in the central overlap volume of Au+Au 
collisions - within UrQMD is about 0.07 GeV/fm 3 at 10.7 
A GeV, 0.1 GeV/fm 3 at 40 A GeV and 0.15 GcV/fm 3 
at 160 A GeV (cf. Ref. j^). These numbers should 
be compared to lattice QCD calculations. Adopting the 
results from Ref. |7j and using T c w 172 MeV (follow- 
ing Ref. we end up with a pressure (in thermal 
and chemical equilibrium) of 0.15 GeV/fm 3 at vanish- 
ing baryon density (jjt B = 0), 0.17 GeV/fm 3 (fi B = 0.21 
GeV), 0.2 GeV/fm 3 (fi B = 0.33 GeV), 0.22 GeV/fm 3 
(/is = 0.42 GeV), and 0.26 GeV/fm 3 (fi B = 0.53 GeV) 
at a temperature T = 1.1T C (~ 190 MeV). These values 
are higher than those from the transport calculations, 
but correspond to a temperature significantly above T c . 
The important finding from the lattice QCD calculations 
is that the pressure drops by an order of magnitude 
when the temperature decreases to T c , which leads to 
~ 0.02-0.03 GeV/fm 3 at T c . These numbers now are 
lower than the transport results for bombarding ener- 
gies > 10.7 AGeV! Consequently, the phase transition 
should be reached even at a lower bombarding energy in 
line with our findings in Fig. 3. 

The authors of Ref. j2^ have proposed to interpret the 
approximately constant slopes above ~ 30A-GeV as 
an indication for a phase transition p9l |. However, to 
support their arguments more experimental data in the 
energy range 20 GeV < y/s < 100 GeV are needed to 
clarify this issue experimentally. 

Summarizing this study, we conclude that according to 
two independent string-hadron transport models, HSD 
and UrQMD, the inverse slope parameters T for K ± 
mesons are practically independent of system size from 



pp up to central Pb+Pb collisions and show a slight in- 
crease with collision energy. The calculated transverse 
mass spectra are in reasonable agreement with the exper- 
imental results for pp collisions and central collisions of 
light nuclei (C+C and Si+Si). The rapid increase of the 
inverse slope parameters of kaons for collisions of heavy 
nuclei (Au+Au) found experimentally in the AGS en- 
ergy range, however, is not reproduced by both models 
(see Fig. [3Jl . Since the pion transverse mass spectra are 
described sufficiently well at all bombarding energies, the 
failure has to be attributed to a lack of pressure especially 
for the strange .K^-mesons. We have argued - based on 
lattice QCD calculations for the pressure at finite tem- 
perature and baryon chemical potential [i B from Ref. 
- that this additional pressure must be generated in the 
early phase of the collision because the strong hadronic 
interactions in the later stages don't produce it. 

The interesting finding of our analysis is, that non- 
hadronic degrees of freedom seem to play a substantial 
role in central Au+Au collisions already at AGS energies 
above ~ 5 A- GeV. Further indications for a low energy 
onset of a phase transition are the missed K + /ir + ratios 
in the AGS and lower SPS regime ^t|- Additionally, the 
excitation function of the elliptic flow t>2 indicates a soft- 
ening of the equation of state at ~ 5 A-GeV in Au+Au 
collisions 30. .31], which signals a change of the pressure, 
too. In spite of the latter observations - all showing a sig- 
nificant change at ~ 5 A-GeV Au+Au collisions - data 
for weakly interacting (on the hadronic level) 0-mesons 
will be necessary to clarify further the question raised in 
the title from the experimental side. 
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